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Abstract 

Our several recent improvements of the current textbook (so called Rayleigh-Schr6ding< 
perturbative representation of bound states are reviewed. They were all inspired by 
an adaptive re-split H(\) = + XH^ of the Hamiltonian. Their feasibility is 
facilitated by the use of nonstandard bases and by the flexibility of normalization of 
the wave functions ijj(X). 



1 Introduction 



Anharmonic oscillators often serve as a methodical laboratory in quantum mechan- 
ics PJ. Schrodinger equation which determines their bound states is, in the time- 
independent case, an elliptic partial differential boundary- value problem, currently 
analysed by the variational or perturbation techniques. In the latter context, one of 
the most transparent illustrations of the implementation and of the efficiency of the 
method is provided by the elementary one- dimensional quartic example 

d 2 

-—^(x) + V(x)ifj(x) = Ei(j(x) (1) 
dx A 

with boundary conditions -0(±oo) = 0, one-parametric family of potentials V(x) = 
x 2 + A x 4 and with the so called Rayleigh-Schrodinger formal-series solution 

E = E<® + A£« + \ 2 EM + ... . 

The example nicely illustrates the main merits of the Rayleigh-Schrodinger approach 
(feasibility, numerical efficiency near A ~ 0) as well as several of its most character- 
istic difficulties. 

In particular, the mere asymptotic and, hence, divergent-series character of the 
anharmonic-oscillator "solutions" (0) divided the physics community into supporters 
and opponents of the perturbative considerations. During the last thirty years of 
development, the former group definitely prevailed. In what follows, an illustration 
of the latter tendency will be presented via a review of some of a few of my own 
recent innovations of the general perturbation formalism. 
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2 Formalisms which do not need the unperturbed 
basis 

One of the most important merits of perturbative considerations lies in their an- 
alytic treatment of certain parameters and in their enormous (though not always 
fully appreciated) formal flexibility. In this sense, there exists an ample variety of 
alternatives to the standard textbook representation of wave functions (0) in the 
(orthonormalized) basis { \n) } of the A = eigenstates of H(0), 

oo 

^ k \x) = J2(x\n)pW . (3) 

n=0 

In this review of my own contributions to this effort I would like to mention the 
representations of wave functions 

• employing the Runge-Kutta-like discretizations of the domain of coordinates, 
both in one and more || spatial dimensions; 

• using the bases with nontrivial orthogonality: Sturmianic [Q, complexified || 

etc; 

• switching to techniques which do not need any scalar product at all and com- 
prise the Taylor-like expansions in the origin || as well as near infinity . 

All these alterations of eq. (|3|) are able to simplify, first of all, the evaluation of the 
matrix elements of the perturbation Hamiltonians H^ k \ 

In this context one should also mention a closely related trick of working with the 
so called "whacko" , order-dependent auxiliary Hamiltonians and bases as employed, 
e.g., in the non-linear context of the so called delta-expansions ||. 
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3 Formalisms with non-diagonal propagators 

In the purely technical sense, the appeal of the textbook perturbative considerations 
may be traced back to their frequent use of the fully diagonalized zero-order models 
H < -°- ) . In practice, an availability of such a type of an approximant is in fact very rare. 
At the same time, people are often seduced by the extremely transparent formulae 
obtainable due to the diagonality of the necessary (so called unperturbed) auxiliary 
propagators R (schematically: R ~ (H^ + const) 

Here, we would like to emphasize, as one of the basic ideas of the present review, 
that the construction of the perturbation series (0) does not necessarily require a 
complete solvability /solution of the Schrodinger equation 

H(X) |V(A)> = E(X) m\)) (4) 

at any particular (usually, zero) value of the auxiliary parameter A = A^ -*. 

There exists, in fact, a balance between the strongly redundant [[| requirements 
of a complete solution of the zero-order problem 

(#(A (0) ) -£(A (0) ))|?KA (0) )> = (5) 

and the related enhancement of economy of the solution of the first-order problem 

(#(A (0) ) - £(A (0) ))|V {1) > + (# (1) - £ (1) )|V>(A (0) )> = (6) 

and of its further higher-order descendants 

(H(X {0) ) - E(X i0) ))\^ k) } = £( fc V(A (0) )> + |r). (7) 

One has to defend a compromise and insist on the obvious observation that the inter- 
related pair (or rather multiplet of the K + 1 linear operator equations) (|) and (0), 
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i.e., the set of the linear operator equations 

A(E®)x = O t A{E^)^ = P- 1 \E^) } k = l,2,...,K (8) 

should be considered as a single, unseparated (and, in general, purely numerical) 
problem of determination of the unknown quantities x, E^ -1 ^ and y^. 

In this sense, one can accept various computational strategies. In a way depending 
on our deeper insight in the structure of the operator A and in its explicit computa- 
tional representation, one could propose the use of some "surviving" simplifications 
offered, e.g., 

• by any form of an incomplete analytic solvability of the zero order equation (El); 
whenever available, such a solvability is still able to simplify the computations 
significantly [|Kj; 



by any algebraization of the Schrodinger eq. (|]); it provides a new origin of 
perturbative constructions as exemplified by the non- variational recipe of ref. 
TI[] or, in the much more numerical context, by the algorithms of ref. ||12|| . 



In all the above-mentioned examples a key role is played by the band-matrix structure 
of the unperturbed Hamiltonian. Also the sparse-matrix structure of its perturba- 
tive "user-friendly" components can make the resulting solution of the whole set 
of equations @ very straightforward and comfortably feasible. Such a construction 
proceeds in fact along the same lines (and with almost the same ease) as in our intro- 
ductory anharmonic example, in spite of the fact that the unperturbed propagators 
themselves remain significantly off-diagonal. 
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4 Working without the model space projectors 

The main price we have to pay for the enhanced flexibility of the above innovated 
algorithms is that we lose an immediate insight in the structure of the unperturbed 
problem itself. 

Several tricks may help to elucidate this structure, at least, indirectly: 

• New attention can be paid to the large but still rigorously terminating cases 
where, i.a., a new, triple-series formalism has been shown manageable ]T5|; 



More attention can be paid to the optimization of our choice of the zero-order 
Hamiltonian operator itself. For this purpose, the above-mentioned Runge- 
Kutta perturbation theory has been combined with an iterative improvement 
of # (0) , e.g., in ref. WA . 



• In an active manner, we may get even much further. Thus, controling the 
mutual methodically relevant relationship between the feasibility of solution 
of the homogeneous and non-homogeneous parts of the Rayleigh-Schrodinger 
hierarchy of equations (|§|) we can introduce also an ad hoc, Hartree-Fock-style 
re-arrangements of many error terms, etc. More technical details about this ef- 
ficient strategy can be found, e.g., in ref. |15[ as well as in the further references 
listed therein. 

One of the main merits of our various modified perturbation prescriptions may be 
seen in their shared enhanced flexibility and in their consequently recurrent character. 
Nevertheless, the related algorithms minimizing the distance between H(X) and 
also exhibit a serious weak point: By their construction, they would all fail in the 
(quasi-) degenerate systems. In this context, the most promising development seems 
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to be represented by the following two different strategies. 

• One can get rid of the model space projectors completely. The first attempt 
in this direction is provided by an inverse-iteration-inspired new scheme of ref. 
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One can employ the shooting-like philosophy as implemented in full detail in 



ref. 17 . 



Technically, a key problem lies in a preservation of feasibility of evaluation of the sep- 
arate Rayleigh-Schrodinger perturbation coefficients. In this setting the relaxation 
of the usual orthogonality of the bases and Hermiticity of the Hamiltonians seems 
to open an access to many new unperturbed H^ '. Simultaneously, a weakening and 
reduction of the traditional exact solvability requirements opens an access to many 
new partially or numerically solvable zero-order approximants with an immediate 
approximative role in realistic models. 



5 Outlook 

The main common purpose of the present proposal of transition to the semi- or 
numerical evaluation of corrections lies in the related possibility of a consequent 
minimization of the absolute magnitude of the perturbation corrections themselves. 
In this sense, the "size" (whatever it means) of the perturbation Hamiltonian is 
efficiently suppressed via a weakening of the technical constraints imposed upon its 
eligible zero-order partners. 

We may summarize that the measurable characteristics of various systems in 
nuclear and atomic physics or quantum chemistry may be fairly well approximated 
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using some of the reviewed non-standard perturbation methods. They reach beyond 
the vicinity of the usual exactly solvable zero order models. This significantly extends 
the domain of applicability of perturbation expansions to many new interesting and 
comparatively complicated systems. 
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